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The  segmented-in-series  solid  oxide  fuel  cell  comprising  fuel  channel,  anode,  cathode  and  electrolyte 
layers  has  been  evaluated  by  developing  a  two-dimensional  model,  in  which  the  equations  have  been 
solved  numerically  through  finite  element  methods.  The  results  indicate  that  the  voltage  of  each  mem¬ 
brane  electrode  assembly  (MEA)  exhibits  a  parabola-like  curve  and  is  higher  than  the  appointed  voltage 
of  unit  cell  (0.7  V).  From  fuel  inlet  to  outlet,  the  voltage  of  each  MEA  deceases  due  to  the  decreasing  local 
H2  concentration.  When  both  the  interconnector  and  electrolyte  gap  lengths  are  fixed,  the  cell  module 
with  5  mm  long  anode  gives  the  maximal  power  density  for  the  SS-SOFC.  Higher  power  densities  can  be 
achieved  through  increasing  the  cathode  thickness. 

©  2009  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Solid  oxide  fuel  cell  (SOFC)  is  one  of  the  most  promising  can¬ 
didates  for  power  generation  because  of  its  high  efficiency,  clean 
generation  of  electric  power  and  silent  work  [  1  -3  ].  The  segmented- 
in-series  SOFC  (SS-SOFC)  has  a  short  current  path  compared  with 
tubular  SOFC  and  the  less  need  for  gas-tight  seals  relative  to  planar 
SOFC  [4-6].  Recently,  SS-SOFC  has  been  implemented  in  a  circular 
tubular  and  a  flat  tubular  configuration.  Several  different  schemat¬ 
ics  of  SS-SOFC  were  proposed  in  different  literatures  [4-10].  A 
fundamental  principle  ought  to  be  clarified  to  avoid  the  short  cir¬ 
cuit  current  inside  the  cell.  Such  as,  the  short  circuit  current  can 
occur  from  the  design  of  Fig.  1  [4]  because  oxygen  anion  can  flow 
through  two  paths  according  to  the  directions  of  two  arrows  shown 
in  Fig.  1. 

As  ohmic  losses  account  for  a  large  part  of  the  overall  losses,  the 
previous  investigations  [4-6]  mainly  focus  on  ohmic  losses  in  SS- 
SOFC.  Besides  the  estimation  of  the  influence  of  ohmic  resistance 
on  cell  performance,  some  other  phenomena  should  also  be  ana¬ 
lyzed  in  detail  for  a  better  understanding  of  this  type  of  cells.  For 
examples,  the  effects  of  both  ohmic  losses  and  fuel  concentration 
distribution  on  the  voltage  of  unit  cell  in  the  SS-SOFC.  An  extensive 
model  including  activation  and  concentration  effects  should  also 
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be  given  for  cell  geometrical  optimization  through  finite  element 
methods  (FEM). 

In  the  present  work,  a  tubular  SS-SOFC  comprising  fuel  chan¬ 
nel,  anode,  cathode  and  electrolyte  layers  will  be  evaluated  by 
developing  a  two-dimensional  (2D)  model.  The  model  covers  the 
electrochemical  reaction,  electrical  and  species  transport,  and  the 
equations  in  this  model  can  be  solved  through  FEM.  The  distribu¬ 
tions  of  electrical  potential  and  fuel  concentration  in  the  cell  will 
be  estimated.  The  mathematical  model  will  be  used  to  carry  out  the 
improvement  of  the  cell  geometry. 

2.  Numerical  model 

Fig.  2  shows  the  schematic  of  SS-SOFC  geometry  in  2D  model. 
The  2D  model  is  axisymmetric,  and  the  dash-dotted  line  is  the  sym¬ 
metric  axis.  H2  humidified  at  room  temperature  is  introduced  as  a 
fuel  at  a  volume  flow  rate  of  30  cm3  min-1.  Oxygen  concentration 
at  the  cathodic  surface  is  assumed  to  be  the  same  as  in  air.  The 
models  of  cell  module  of  the  different  dimension  are  established 
to  investigate  the  effect  of  the  cell  structure  on  cell  performance. 
Fig.  2  shows  a  schematic  of  the  repetitive  unit  of  the  SS-SOFC  as 
well.  It  is  assumed  that  interconnector  and  electrolyte  gap  lengths 
are  fixed  at  all  the  process  of  calculations.  The  interconnector  and 
electrolyte  gap  lengths  are  1 .5  and  2  mm,  respectively.  The  mate¬ 
rial  properties  are  listed  in  Table  1 .  The  SS-SOFC  model  chooses  the 
electrolyte  thickness  of  10  pum  and  the  anode  thickness  of  20  pim. 
The  cathodic  thickness  and  the  inner  radius  of  the  anode  chan¬ 
nel  are  changed.  Because  of  the  cathodic  conductivity  being  lower 
than  the  anodic  conductivity,  it  is  considered  to  change  the  cathodic 
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Nomenclature 

c  concentration  (mol  rrr3 ) 

Dy  diffusion  coefficient  (m2  s-1 ) 

F  Faraday’s  constant  (96,485  C  mol-1 ) 

i  local  current  density  (A  m-2 ) 

ne  electron  transferred  per  reacting 

P  pressure  (Pa) 

Rg  gas  constant  (8.314 J  mol-1  K-1 ) 

T  temperature  (I<) 

u  velocity  (m  s-1 ) 

Greek  letters 

P  transfer  coefficient 

p  density  (kg  rrr3) 

a  conductivity  (S  nrr 1 ) 

0  potential  (V) 

Subscripts  and  superscripts 
a  anode 

c  cathode 

ele  electrolyte 

rev  reversible 


thickness  rather  than  the  anodic  in  this  paper.  The  thickness  of  the 
support  substrate  of  the  SS-SOFC  is  1  mm.  The  total  length  of  cell 
module  is  about  40  mm,  which  depends  on  both  the  anode  length 
and  number  of  unit  cell.  The  cell  dimensions  of  the  different  anode 
length  are  shown  in  Table  2. 

In  this  steady  state  model,  the  reactant  gas  mixtures  are  approx¬ 
imated  as  ideal  gas.  The  electrochemical  reactions  are  considered 


Fig.  1.  Schematic  drawing  of  the  SS-SOFC  geometry  in  Ref.  [4]. 


Table  1 

Input  parameters  to  model  [2,3,11]. 


Descriptions 

Symbol 

Value 

Anode  thickness  (|xm) 

20 

Electrolyte  thickness  (p,m) 

10 

Anode  conductivity  (S  m-1 ) 

95  x  106  T_1  exp(-1150/T) 

Cathode  conductivity 

(Sm-1) 

crc 

42  x  106  T-1  exp  (-1200/T) 

Electrolyte  conductivity 
(Sm-1) 

CTeie 

3.34  xlO4  exp  (-10,300/7) 

Interconnector 
conductivity  (S  itt1  ) 

Vic 

40 

Transfer  coefficient 

0.5 

Activation  energy  of  anode 
(J  mol-1 ) 

Fact  ,a 

140,000 

Electron  transferred  per 
reacting 

ne 

1 

Hydrogen  diffusion 
coefficient  (m2  s_1) 

Du2,a 

3.14e-5 

H20  diffusion  coefficient 

(m2  s-1) 

Dn2o,a 

1.39e-5 

Gas  diffusion  coefficient  in 
anodic  channel  (m2  s_1 ) 

Da 

8.506e-4 

Oxygen  diffusion 
coefficient  (m2  s_1) 

Do2,c 

7.588e-6 

e 


Fig.  2.  Schematic  drawing  and  repetitive  unit  of  a  section  of  the  SS-SOFC  geometry. 
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Table  2 

Cell  dimensions  of  different  anode  length. 


Anode  length  (mm) 

Unit  cell  length  (mm) 

Unit  cell  number 

Active  length  of  cell  module  (mm) 

Total  length  of  cell  module  (mm) 

3 

1.5 

7 

10.5 

42 

4 

2.5 

6 

15 

41.5 

5 

3.5 

5 

17.5 

39 

7 

5.5 

4 

22 

38.5 

10 

8.5 

3 

25.5 

37 

18 

16.5 

2 

33 

40 

to  take  place  only  at  the  electrode/electrolyte  interface.  The  model 
is  assumed  to  be  isothermal  and  the  cell  module  runs  at  800  °C. 
It  is  assumed  that  the  average  voltage  of  each  unit  cell  is  0.7  V 
for  the  operation  of  SS-SOFC.  For  example,  when  there  are  7  unit 
cells,  the  cell  module  voltage  is  7  x  0.7  =  4.9  V.  The  model  makes 
use  of  the  differential  equations  which  are  integrated  in  the  2D 
domains  through  commercial  software,  COMSOL  MULTIPHYSICS®. 
Input  parameters  of  the  model  are  shown  in  Table  1. 

The  Butler-Volmer  kinetic  equations  are  generally  adopted  for 
expressing  the  local  current  density  distribution: 


j  _  jH,  (Ph2  )C 

la  ~  la,0 


(P  H20 


,0.5 


exp 


—Eact,a 

~R^r 


x  <  exp 


P  ne  •  F  •  (0rev  10a  0ele  I )  A 


Rg-T 


-  exp 


-(W) 


ne  F  •  (0rev  —  10a  —  0e/el) 


Rg-T 


(1) 


Fig.  4.  The  voltage  (V)  of  each  MEA  for  the  different  anode  length  with  the  support 
substrate  of  0.8  mm  inner  radius  and  the  cathode  of  50  pirn  thickness. 


k  =  %  ■ 


(Po2/Pq2)1/4 


1  +(Po2/pS2)V2 


x  <J  exp  — l0c-0elelA 


-  exp 


-(1-0) 


Tle  F  ■  (0rev  —  10c  —  0e/e  I ) 


Po,  (Pa)  =  4.9  X  103  exp  (  - 


2  x  105 


In  Eqs.  (1)  and  (2a),  the  anodic  exchange  current  density, 
is  3.4  x  1010  Am-2  for  a  Ni-YSZ  anode  [12].  The  cathodic  exchange 
current  density,  2  is  4000  Am-2  for  a  LSM-YSZ  cathode  [3].  The 
values  of  ijjh  and  ij? 2  have  been  used  in  the  FEM  model  of  SOFC 
[3],  where  the  results  showed  that  the  model  agrees  well  with  the 
experimental  data.  0rev  can  be  obtained  from  Nernst  equation.  Elec¬ 
tron  and  ion  charge  conservation  at  the  electrode  and  electrolyte  is 
(2a)  expressed  by  Eq.  (3): 

-V(cr  •  V0)  =  0  (3) 

(2b)  The  steady  state  diffusion  and  convection  equation,  Eq.  (4)  is 

solved  to  obtain  the  species  distribution  of  reactants  and  products 


Potential  (V)  Current  density  (A  nr2) 


Fig.  3.  (a)  Potential  (V)  distribution  of  the  SS-SOFC.  (b)  Streamline  and  distribution  plots  of  current  density  (Am-2)  from  the  cathode  via  the  interconnector  to  anode. 
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Fig.  5.  The  voltage  (V)  of  the  second  MEA,  the  potentials  of  its  cathode  and  anode  for 
18  mm  anode  length  with  support  substrate  of  0.8  mm  inner  radius  and  the  cathode 
of  50  |jim  thickness. 


in  the  gas  channel  and  the  electrodes: 

V(-Dy  •  Vq  +  Cj  •  u)  =  0  (4) 

Anode/electrolyte  interface  boundary: 

(-^H2,aVcH2  +  Ch2u)  •  n  =  (4a) 

(-^H2,aVCH2o  +  Ch2oU)  *  H  ^  (4b) 

Anode  channel  inlet: 

Ch2  =  Ch2,o  (4c) 

Ch2o  =  Ch2o,o  (4d) 


(a)  (b)  (c)  (d)  (e)  (f) 

Fig.  6.  Surface  plots  of  H2  concentration  (molm-3)  distribution  for  the  different 
anode  length  of  the  SS-SOFC  (anodic  length  from  (a)  to  (f)  is  3, 4,  5, 7, 10  and  18  mm, 
respectively). 


Fig.  7.  H2  concentrations  (molm-3)  along  the  anode/substrate  interface  and  fuel 
channel  center  in  the  SS-SOFC  of  the  different  anode  lengths. 


Cathode/electrolyte  interface  boundary: 

(-£>o2,cVco2  +Co2u)-n  =  (4e) 

02  concentration  at  cathode/air  interface: 

Co2  =  Co2,o  (4f) 

The  governing  equations  are  solved  with  the  commercial  soft¬ 
ware,  COMSOL  MULTIPHSICS.  In  this  article,  we  choose  to  create  a 
free  mesh  consisting  of  triangular  elements  for  a  2D  geometry  in 
COMSOL  Multiphysics.  In  these  calculations,  mesh  consists  of  about 
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Fig.  8.  Effects  of  different  geometries  on  power  density  of  the  SS-SOFC. 


140,000  elements.  The  convergence  criterion  is  relative  tolerance 
ofle-6. 

3.  Results  and  discussion 

3.1.  Transport  phenomena  analysis  in  SS-SOFC 

The  potential  distribution  of  the  cell  module  is  shown  in  Fig.  3. 
The  potential  is  highest  on  the  right-hand  side  where  the  cathode 
connects  with  the  current  collector.  The  lowest  potential  is  on  the 
left-hand  side  where  the  anode  connects  with  the  current  collector. 
As  far  as  the  cell  module  is  concerned,  the  potential  deceases  from 
right  to  left.  For  the  change  of  the  local  potential  in  Fig.  3,  it  is  shown 
that  the  potential  deceases  from  left  to  right  inside  the  rectangle 
and  the  arrow  shows  the  current  direction.  At  the  connecting  region 
of  two  unit  cells,  the  potential  decease  from  the  cathode  via  the 
interconnector  to  anode. 

The  distribution  of  the  current  density  is  not  uniform  in  the  unit 
cell.  To  comprehend  the  flow  of  the  current  in  the  SS-SOFC,  the 
streamline  and  distribution  plots  of  current  density  (Am-2)  from 
the  cathode  via  the  interconnector  to  anode  are  shown  in  Fig.  3.  The 
part  enclosed  by  the  rectangle  in  Fig.  3(a)  is  scaled  up  to  Fig.  3(b). 
The  streamlines  express  current-flow  directions  and  paths.  The  cur¬ 
rent  from  the  cathode  to  anode  flows  through  interconnector  with 
the  shortest  path,  which  is  because  the  electrical  conductivity  of 
the  interconnector  is  far  lower  than  those  of  the  electrodes.  It  can 
be  seen  that  the  highest  current  density  occurs  at  this  place  where 
the  path  of  the  current  is  shortest  from  the  cathode  to  the  anode. 
The  current  flowing  from  cathode  to  anode  is  mainly  via  this  area. 
The  phenomena  show  that  the  most  current-flow  from  cathode  to 
anode  crosses  the  interconnector  from  a  very  small  area. 

Fig.  4  shows  the  voltage  of  each  membrane  electrode  assembly 
(MEA)  with  the  support  substrate  of  0.8  mm  inner  radius  and  the 


cathode  of  50  fxm  thickness.  The  anode  length  is  changed  as  3, 4, 5, 
7,  10  and  18  mm.  The  voltage  of  each  MEA  is  the  potential  differ¬ 
ence  between  the  cathode  and  the  anode.  The  voltage  of  each  MEA 
is  higher  than  0.7  V  because  of  the  ohmic  loss  from  interconnec¬ 
tor.  Fig.  4  depicts  that  the  voltage  of  each  MEA  decreases  from  fuel 
inlet  to  outlet  because  the  fuel  concentration  decreases.  The  cur¬ 
rent  flowing  through  each  MEA  in  the  SS-SOFC  is  same.  Flowever, 
the  voltage  is  different  at  each  MEA  because  the  Nernst  voltage  of 
each  unit  cell  is  different.  When  the  longer  SS-SOFC  is  designed, 
the  voltage  is  much  lower  at  the  last  MEA  than  the  front  MEA.  The 
anode  of  the  last  unit  cell  is  at  the  relatively  high  oxygen  pressure. 

Fig.  4  also  depicts  that  the  voltage  of  each  MEA  exhibits  a 
parabola-like  curve.  To  understand  this  phenomenon,  the  voltage 
of  the  second  MEA  with  the  anode  of  18  mm  length  in  Fig.  4  is  cho¬ 
sen.  The  chosen  voltage  curve  is  enclosed  with  the  black  half  circle 
in  Fig.  4.  The  potentials  of  its  cathode  and  anode  are  also  showed 
in  Fig.  5.  It  indicates  that  the  potentials  of  both  cathode  and  anode 
increase  from  fuel  inlet  to  outlet  and  express  a  protruding  and  con¬ 
cave  curve,  respectively.  The  voltage  of  each  MEA  is  the  potential 
difference  between  cathode  and  anode,  so  the  voltage  curve  of  each 
MEA  is  parabola. 

The  distributions  of  hydrogen  concentration  for  the  SS-SOFC 
of  the  different  anodic  length  are  shown  in  Fig.  6.  Since  the  con¬ 
sumption  of  hydrogen  occurs  at  the  anode/electrolyte  interface, 
the  hydrogen  concentration  decreases  along  the  fuel-flow  direc¬ 
tion.  Fig.  6  shows  that  hydrogen  concentration  is  low  close  to  each 
anode.  Fig.  7  shows  the  variation  of  H2  concentrations  (molm-3) 
along  both  the  anode/substrate  interface  and  fuel  channel  cen¬ 
ter  of  the  SS-SOFC  for  the  different  anode  lengths.  The  legend 
“3  mm-central”  means  the  variation  of  FI2  concentrations  along 
fuel  channel  center  of  the  SS-SOFC  for  the  anode  length  of  3  mm. 
The  legend  “4  mm-interface”  means  the  variation  of  H2  concen¬ 
trations  along  the  anode/substrate  interface  of  the  SS-SOFC  for  the 
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anode  length  of  4  mm.  The  H2  concentrations  of  all  cell  modules  are 
the  same  at  inlet  and  decreases  along  the  fuel-flow  direction.  The 
local  H2  concentration  at  the  anode/electrolyte  interface,  where 
H2  is  consumed,  is  locally  minimum.  The  H2  concentrations  show 
a  sawtooth-like  behavior,  with  the  period  set  primarily  by  the 
anode  length.  A  less-frequency  oscillation  can  be  seen  when  the 
anode  length  increases.  The  variation  of  H2  concentrations  along 
fuel  channel  center  is  similar  to  the  H2  concentrations  at  the  active 
interface,  but  its  variation  is  quite  weak.  Comparing  Fig.  4  with 
Fig.  7,  it  is  clearly  seen  that  the  variation  of  the  local  FI2  concen¬ 
tration  causes  the  decrease  of  the  voltage  of  each  MEA  from  fuel 
inlet  to  outlet. 

3.2.  Improvement  of  cell  geometry 

The  mathematical  model  is  used  to  carry  out  the  improvement 
of  cell  geometry.  The  effects  of  cell  geometry  on  the  cell  perfor¬ 
mance  are  considered  for  the  improvement  of  the  SS-SOFC.  The 
operating  temperatures  of  all  cells  are  800  °C  in  the  calculations. 
Fig.  8  displays  the  effects  of  the  anode  lengths  on  the  power  den¬ 
sities  of  the  SS-SOFC.  In  this  figure,  power  densities  are  calculated 
using  power  divided  by  the  length  of  the  cell  module.  The  power 
densities  decrease  at  longer  unit  cell  lengths  due  to  high  resistance 
loss,  whereas  they  decrease  at  shorter  unit  cell  lengths  due  to  less 
unit  cell  active  area.  The  maximal  power  density  is  a  result  of  the 
influences  of  both  electrode  resistance  and  unit  cell  active  area  with 
the  change  of  the  anode  length.  The  optimal  length  of  the  anode  is 
5  mm,  namely,  the  active  length  of  unit  cell  is  3.5  mm. 

Besides  the  anode  length,  the  other  key  factors  affecting  cell 
performance  are  the  cathodic  thickness  and  the  inner  radius  of  the 
fuel  channel.  Fig.  8  shows  the  effects  of  both  the  cathode  thickness 
and  the  inner  radius  of  the  fuel  channel  on  the  power  densities  for 
the  different  anodic  lengths.  An  increase  of  the  cathode  thickness 
leads  to  a  reduction  of  ohmic  losses  of  the  SS-SOFC,  and  causes 
a  significant  increase  of  the  cell  performance.  The  increment  of 
power  density  is  very  little  when  the  cathode  thickness  is  larger 
than  75  p,m,  because  the  losses  of  the  concentration  diffusion  play 
an  important  role  for  such  a  thick  cathode.  An  increase  of  the  inner 
radius  of  fuel  channel  can  lead  to  an  increase  of  cell  module  active 
area  in  the  SS-SOFC,  and  thus  raise  the  power  densities. 


4.  Conclusions 

The  2D  mathematical  models  were  developed  to  predict  the 
electrical  and  species  transport  behavior  in  the  SS-SOFCs  of  differ¬ 
ent  geometries.  The  results  indicate  that  the  voltage  of  each  MEA 
exhibits  a  parabola-like  curve  and  is  higher  than  the  averaged  volt¬ 
age  of  unit  cell  (0.7  V).  The  voltage  of  each  MEA  deceases  in  the 
fuel-flow  direction  due  to  the  decrease  of  the  local  FI2  concentra¬ 
tion.  When  the  electrolyte  and  interconnector  gaps  are  fixed,  the 
cell  module  with  the  anode  length  of  5  mm  gives  the  maximum 
of  power  density  for  the  SS-SOFC.  Fligher  power  densities  can  be 
obtained  through  increasing  the  cathode  thickness.  These  calcu¬ 
lations  can  be  used  to  design  and  optimize  SS-SOFC  for  the  given 
materials  and  operating  condition. 
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